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Hofstadter butterfly for a finite correlated system 
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We investigate a finite two-dimensional system in the presence of external magnetic field. We 
discuss how the energy spectrum depends on the system size, boundary conditions and Coulomb 
repulsion. On one hand, using these results we present the field dependence of the transport prop- 
erties of a nanosystem. In particular, we demonstrate that these properties depend on whether the 
system consists of even or odd number of sites. On the other hand, on the basis of exact results 
obtained for a finite system we investigate whether the Hofstadter butterfly is robust against strong 
electronic correlations. We show that for sufficiently strong Coulomb repulsion the Hubbard gap 
decreases when the magnetic field increases. 



I. INTRODUCTION 

The problem of electrons moving in a periodic potential 
under the influence of an external magnetic field has been 
investigated since the beginning of quantum mechanics. 
Despite the seeming simplicity of the problem, many its 
aspects still remain unresolved. Even in the absence of 
electronic correlations solutions are known only in limit- 
ing cases. In particular, two dimensional (2D) electron 
gas under the influence of a periodic potential and a per- 
pendicular magnetic field can be described in two limits, 
one of a weak and the other of a strong periodic poten- 
tial. In the former case the applied magnetic field is the 
main factor that determines the behavior of electrons. As 
a result the electronic wave functions are Landau level- 
like, with the degeneracy lifted by the periodic poten- 
tial. If the potential is modulated in one dimension, the 
width of the resulting "Landau bands" oscillates with the 
magnetic field as a consequence of commensurability be- 
tween the cyclotron diameter and the period of potential 
modulation. It leads to oscillations in magnetoresistance, 
known as the Weiss oscillations. 1 If the potential is modu- 
lated in two dimensions, "minigaps" open in the "Landau 
bands" , and the energy spectrum plotted versus the ap- 
plied field composes the famous Hofstadter butterflyi^ It 
is interesting, that the same spectrum occurs in a comple- 
mentary limit, when the lattice potential is very strong, 
and the electronic wave functions are Bloch-like, modi- 
fied by the magnetic fields 

The simplest model for the case, when an applied 
field and a lattice potential are present simultaneously, 
is commonly referred to as the Hofstadter or Azbel- 
Hofstadter modeli^ The corresponding Hamiltonian de- 
scribes electrons on a two-dimensional square lattice with 
nearest-neighbor hopping in a perpendicular uniform 
magnetic field. The Schrodinger equation takes the form 
of a one-dimensional difference equation, known as the 
Harper equation (or the almost Mathieu equation) 
It is also a model for a one-dimensional electronic sys- 
tem in two incommensurate periodic potentials. The 
Harper equation also has links to many other areas of 
interest, e.g., the quantum Hall effect; 7 - quasicrystals, 
localization-delocalization phenomena;^ the noncom- 



mutative geometry^ the renormalization group ; ' the 
theory of fractals, the number theory, and the functional 
analysis^ It is also useful in determining the upper crit- 
ical fiel d 14 ' 15 ! 16 ' 17 and the pseudogap closing fields in 
high-temperature superconductors . 

The unusual structure of the Hofstadter butterfly is 
a characteristic feature of 2D systems. Similarly to the 
case of Landau levels the dimensionality of the system is 
of crucial importance for the Hofstadter butterfly. Move- 
ment of electrons along the external magnetic field would 
be responsible for broadening of the Hofstadter bands. 
They eventually may overlap and, in this way, smear 
out the original fractal structure of the energy spec- 
trum. On the other hand it is known, that electronic 
properties of low dimensional systems may be completely 
changed by the presence of electronic correlations. In 
particular, it is well known for ID systems that pertur- 
bation theory breaks down and arbitrarily weak on-site 
Coulomb repulsion qualitatively changes the whole exci- 
tation spectrum from the Fermi to Luttinger liquid type. 
The role of Coulomb interaction in the case of 2D lat- 
tice has intensively been investigated in connection with 
high temperature superconductors. Although a complete 
description of the correlated 2D system is still missing, 
it became obvious that mean-field approaches are in- 
valid even for moderate values of Coulomb repulsion. 
Therefore, a question arises, whether the fine structure 
of the Hofstadter butterfly is robust against the presence 
of electronic correlations. This problem has previously 
been investigated on a mean-field leveL 19 ' 20 In particu- 
lar, the analysis presented in Ref. [2(3 suggests that also 
in the presence of electronic correlations the energy levels 
should form the Hofstadter butterfly with additional en- 
ergy gap in the middle of the energy spectrum. However, 
the mentioned above limitations of the mean-field results 
clearly show, that these results do not represent a con- 
clusive solution. In this paper we address this problem 
with the help of a method that is particularly suited for 
investigations of the low dimensional correlated systems, 
i.e, exact diagonalization of finite clusters. Although in 
this method short-range electronic correlations are ex- 
actly taken into account, finite size effects will seriously 
affect the energy spectrum. These modifications can be of 
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special importance, when the system is under influence 
of magnetic field i 21 ' 22 In order to separate the correla- 
tions and size-induced effects we start our investigations 
with a finite uncorrelated system and discuss both fixed 
(fbc) and periodic boundary conditions (pbc). Apart 
from the discussion of the Hofstadter butterfly these re- 
sults may be applicable to investigations of nanosystems 
in the presence of magnetic field, where fixed boundary 
conditions are more appropriate than the periodic ones. 



II. FINITE-SIZE EFFECTS 

For the sake of completeness, we start with a brief 
derivation of the Harper equation for the gauge A = 
B (—ay, (1 — a)x, 0). Here, the parameter a e (0, 1) al- 
lows one to distinguish between the Landau (a = 0) and 
symmetric gauges (a = 1/2). The 2D square lattice in 
the presence of external, perpendicular magnetic field can 
be described by the tight-binding Hamiltonian: 

x,y,a 

-t r, _l_ J-4>y'-cJ r , . \ 



e -i0„a; c t 



(1) 



where ct , creates an electron with spin a at the 
site (x,y) and t is the nearest-neighbor hopping inte- 
gral in the absence of magnetic field. <f> x = 2iraa, 
(f> y = 2ira(l — a), a = $/3>0; where $ is the magnetic 
flux through the lattice cell and <i>o is the flux quantum. 
In order to determine cigcnfunctions of this Hamiltonian, 
formally one should solve a 2D eigenproblem. However, 
in the case of the Landau gauge (a = 0), the hopping 
integrals in Eq. ([T|) depend solely on the x coordinate. 
Because of the translational invariance along the y axis, 
eigenf unctions exhibit a plane-wave behavior in this di- 
rection [exp(iky)]. This argumentation can be extended 
to a more general gauge. Such an extension requires an 
appropriate shift of the momentum k. For a N x x N y 
lattice with pbc, we introduce fermionic operators c k 
defined by 



"(x,y)<r 




-iy(k+4> x x) t 



^x,k.cr ' 



(2) 



The Hamiltonian (1) can then be written: 



H 



x.k.a 



2cos(fc — 27rax)c x 



h\<T 



(3) 



The resulting Hamiltonian is diagonal in the quantum 
numbers k and tridiagonal in the x-coordinates. It means 
that the applied ID transformation to the momentum 
space, allows one to reduce the original 2D eigenprob- 
lem [Eq. |(TJ)] to a ID one. In the following we study 
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FIG. 1: (Color online). Density of states, p(E), as a func- 
tion of the magnetic field (a) calculated with fixed (left col- 
umn) and periodic (right column) boundary conditions [see 
Eq. @]. The brighter colors correspond to larger values of 
the density of states. The first, second and third row shows 
results obtained for a 6 x 6, 12x12 and 30 x 30 cluster, respec- 
tively. The broadening of the one-particle levels r\ — 0.05t 
has been assumed. In the last row we magnify the regions 
marked above. In order to demonstrate the presence of the 
edge states, the broadening has been reduced to r\ — 0.005t. 



whether such a reduction is possible also for a finite 
system with fixed boundary conditions. The relevant 
wave-function in this case is no longer the plane-wave, 
since it must vanish at the system edges. It is easy 
to check, that the Hamiltonian of a ID chain with fbc 
can be diagonalized with the help of the transformation 
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c y.a ~ Ylk sin(fcy)c|, a . In the case of a 2D lattice in the 
presence of the external magnetic field, an appropriate 
shift of the wave-vectors allows one to cancel the Peierls 
phase factors for the hopping along the y-axis: 

where the wave-vectors 

7T 2-7T NyTT 

~ N y + 1 ' a^TT' '"' a^TT' 

In the above equation we have assumed a = 1. General- 
ization to arbitrary value of a is straightforward. 
The orthogonality relation 

^2 sm(ky)sin(py) = Ny ^ 1 S kp 

y 

allows one to carry out the inverse transformation. Then, 
the Hamiltonian takes the form 

_^ /n,-i w » 2 — , 

H =tJ2( EE^iE si *)™(py) x 

k.a \ x = l y=l « p 

s=l / 

Only the term responsible for the hopping along the y- 
axis is diagonal in k, whereas the remaining hopping term 
is generally not. Therefore, contrary to the case of infinite 
lattice with pbc, the Hamiltonian cannot be reduced to 
a form that is diagonal in wave-vectors and tridiagonal 
in real space coordinates. The only exceptions occur for 
a = and a = 1/2. The first case is trivial. In the latter 
case (a = 1/2), an additional transformation 

Cx,k,a -> ^C X)fc)CT (1 + (-If) - ~C x ^-k, a (1 - (-ID , 

leads to the Hamiltonian in the form given by Eq. |[3]). 
However, since k G (— 7r,7r) for pbc and k € (0, 7r) for 
fbc, also in this case the energy spectrum depends on 
the boundary conditions. Consequently, degeneracy of 
the energy levels is lower for fbc than for pbc. An addi- 
tional difference, that immediately follows from the an- 
alytical calculations, is related to the density of states 
p(e) for a = 1/2. In the case of infinite system with 
pbc p(0) = 0, whereas for a finite lattice with fbc there 
exists an eigenvalue e = 0, provided the system consists 
of odd number of sites. The proof of this statement is 
straightforward and, therefore, we omit the details. The 
difference between systems with even and odd number 
of sites will be discussed in more details in connection 
with the transport properties. In Fig. 1 we compare the 
density of states obtained for various boundary condi- 
tions and lattice sizes. It has been calculated using the 
standard formula: 

p(E) = V3 , (4) 



where e n is the energy of the n-th one-particle eigenstate. 

In the presence of magnetic field the translation group 
of the lattice does not represent a symmetry group of the 
Hamiltonian and one can discuss periodicity only with 
respect to the magnetic translation group i 23 i 24 Conse- 
quently, for a finite system one can apply the pbc only 
for some particular values of a, which are determined 
by the system size. On the other hand, for fbc the en- 
ergy levels can be calculated for arbitrary magnetic field, 
as it was demonstrated in Ref. I22I However, in order 
to compare directly results obtained for fbc and pbc, in 
both the cases we have used only these values of mag- 
netic field, which are allowed for pbc. In the case of 
fbc there exist edge states, which are responsible for ad- 
ditional levels inside the energy gaps in the Hofstadter 
butterfly^ They are clearly visible in Fig. 1 in the den- 
sity of states obtained for a 6 x 6 system. The relative 
contribution of these states to p(e) decreases with the 
system size. Therefore, they are much less visible in the 
density of states obtained for larger systems and become 
unimportant for an infinite lattice. However, for a finite 
system, they may qualitatively change the well known 
structure of the Hofstadter butterfly, obtained from the 
Harper equation. 

A. Transport properties 

The above discussion can directly be applied to the in- 
vestigations of transport properties of nanosystems. The 
field-induced modifications of the energy spectra lead to 
strong changes of the transport current, at least in the 
low voltage regime. We use the formalism of nonequi- 
librium Green functions to analyze these effects in a 
nanosystem coupled to leads. The coupling is described 
by: 

ifnano-el = ^ (9U.,x,y,adka, a C (x,y)cr + H.C.^ , (5) 
k,a:,y,a,er 

where di „ „ creates an electron with momentum k and 
spin a in the electrode a. Here, a € {L,R} indicates 
the left or right electrode. We assume a simple model 
in which the leads are described by a two-dimensional 
(2D) lattice gas and the hopping between the leads and 
nanosystem is possible only perpendicularly to the edge 
of the nanosystem. gk,x,j,,a is nonzero only for sites (x, y) 
which are located at the edge a of the nanosystem. The 
details of calculations can be found in Refs. [26| and [2?]. 
Figs. 2 and 3 show the field dependence of the transport 
current for various size of the nanosystem. 

Two main features arise from the presented results: (i) 
strong magnetoresistance in a weak field regime and (ii) 
even-odd parity effect. Both these features occur for a 
low voltage only. The first effect is shown in Fig. 2. It 
can easily be understood on the basis of the field depen- 
dence of the one-particle energies that are close to the 
Fermi energy, as presented in the inset in Fig. 2. One 
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FIG. 2: (Color online) Field dependence of the transport cur- 
rent through the 5 x 5 and 6x6 nanosystems coupled to macro- 
scopic leads for the applied voltage eV = O.Olt. Iq = 2et/h, 
where e is the elementary charge and h is the Planck con- 
stant. See Fig. 5 in Ref. [2y for parameters and the details of 
the coupling between the nanosystem and the electrodes. The 
inset shows the field dependence of the one-particle eigenen- 
ergies located in the vicinity of E = obtained for an isolated 
nanosystem with fbc. 



can see, that for low voltage the number of states which 
participate in the transport decreases when the magnetic 
field increases. This effect is due to a field-induced split- 
ting of a strongly degenerated level at zero energy. This 
degeneracy, in turn, is a remnant of the van Hove singu- 
larity, that is a typical feature of an infinite 2D lattice. 
There is, however, a significant difference between sys- 
tems with even and odd number of lattice sites. In the 
former case, strong magnetic field can completely remove 
states from the vicinity of the Fermi energy, what results 
in vanishing of the current. On the other hand, if there is 
an odd number of lattice sites, one state is always located 
at zero energy. As a result a finite conductivity occurs 
for arbitrary magnetic field. This parity effect originates 
from the facts that the energy spectrum is symmetric 
with respect to the zero energy and the number of en- 
ergy levels is equal to the number of lattice sites. The 
difference between systems consisting of even and odd 
number of sites is pronounced for a = 1/2, what can be 
seen in Fig. 3. The even-odd parity effect is a well known 
feature of persistent currents in mesoscopic rings, where 
it occurs due to its nontrivial first homotopy group^ 
Here, we have demonstrated that a similar effect may 
occur also in a nanosystem with a trivial topology: in an 
isolated nanosystem as well as in a system coupled to a 
macroscopic leads. In the first case it shows up in the 
energy spectrum, whereas in the latter case it is visible 
also in the transport properties. An additional similarity 
between rings and the systems under investigation con- 
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FIG. 3: (Color online) Upper panel shows the field depen- 
dence of the transport current through the 5x5 and 6x6 
nanosystems for the applied voltage eV = O.Olt. Lower panel 
shows results obtained for 9x9 and 10 x 10 nanosystems with 
eV = 0.002t. The remaining parameters are the same as in 
Fig. 2. 



cerns the fact, that the parity effect occurs only in small 
systems and disappears in the thermodynamic limit. 



III. THE ROLE OF CORRELATIONS 

After establishing the role of boundary conditions and 
its significance for the properties of nanosystems, we 
switch to the main question, whether the Hofstadter en- 
ergy spectrum is robust against the presence of strong 
electronic correlations. Both the external magnetic field 
and electronic correlations give rise to opening of the en- 
ergy gaps in the density of states. One could expect that 
these mechanism should independently contribute to the 
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FIG. 4: Density of states obtained for a 2 x 3 cluster with fbc for various values of the Coulomb repulsion and the magnetic 
field. 



opening of these gaps. In the following, we demonstrate 
that this intuitive statement is wrong and, actually, ex- 
ternal magnetic field reduces the Hubbard energy gap. 
In order to investigate this problem, we consider the 2D 
Hubbard model in the presence of magnetic field: 



HKubb = Hki- 



(6) 



where the kinetic term is given by Eq. |T]) and the elec- 



tron number operator nr Xj 



Z (x,y)cr C ( x ,y) 



In the fol- 



lowing we investigate the half-filled case, when the Hub- 
bard gap opens at the Fermi level, driving the system 
from a metallic state to an insulating one. 

This Hamiltonian has exactly been diagonalized 
mostly by means of the Lanczos method. It is one of 
the most effective computational tools for searching for 
the ground state and some low laying excited states of a 
finite system. We start with a system sufficiently small 
to allow one to determine the whole energy spectrum. 
In the case of finite system calculations pbc can be ap- 
plied only for specific values of the magnetic field, which 
depend on the cluster size. Therefore, in the follow- 
ing we use the fbc. Fig. 4 shows how the density of 
states depends on the magnetic field and the magnitude 
of the Coulomb repulsion. One can see, that for weak 
to moderate Coulomb repulsion its influence on the den- 
sity of states strongly depends on the applied magnetic 



field. In the absence of magnetic field (a = 0) even rel- 
atively weak electronic correlations (U — It") strongly 
modify the energy spectrum. On the other hand, for 
strong magnetic field the density of states becomes ro- 
bust against the Coulomb correlations. In particular, for 
a = 1/2 densities of states obtained for U = and U = 2t 
hardly differ from each other. This result already sug- 
gests that Coulomb interaction modifies various parts of 
the Hofstadter butterfly in a different way. It remains 
in contradiction to the results obtained in the mean-field 
analysis. 20 In the latter approach the Coulomb repulsion 
opens an almost field independent gap in the middle of 
the Hofstadter butterfly. Additionally, the one-particle 
energies remain in the range (—At, At), at least for mod- 
erate U, i.e., the electronic correlations do not change 
the spectrum width. It has also been reported in Ref. [20 
that for stronger Coulomb repulsion the only modifica- 
tions of the Hofstadter butterfly concern the larger band 
gaps and narrower band widths. Contrary to the mean- 
field results, exact diagonalization method indicates that 
electronic correlations lead to a substantial modification 
of the density of states. First, there are excitation beyond 
the range (— At, At), When U increases the one-particle 
excitations are replaced by the collective ones with a fi- 
nite live-time. Instead of the <5— like peaks in the density 
of states we obtained much larger number of wider peaks, 
which eventually may overlap. Therefore, we expect that 
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FIG. 5: (Color online) Magnitude of the reduced Hubbard 
gap (A*(U, a)) as a function of U for various magnetic fields. 
These results have been obtained for a 12-site cluster with 
fbc. 

electronic correlations smear out the fine structure of the 
Hofstadter butterfly. 

For the application purposes, the most important prop- 
erty concerns the density of states in the vicinity of Fermi 
level. Therefore, we restrict the following study only 
to the ground state and the lowest excited states. Es- 
pecially, we analyze field dependence of the Hubbard 
gap A(Z7, a) that opens in the middle of the density of 
states. In Fig 5. we present a reduced gap A*(U, a) — 
A(U,a) — A(U,Q), whereas A(U,0) is presented in the 
inset in this figure. 

As one could expected, A increases with U for arbi- 
trary magnetic field. However, the negative slope of A* 
as a function of U clearly indicates that this increase is 
smaller when the magnetic field is switched on. Finally 
for sufficiently strong Coulomb repulsion, the Hubbard 
gap seems to be a monotonically decreasing function of 
magnetic field for a G (0,1/2). Since the Hamiltonian 
is invariant under the transformation a — > 1 — a the 
gap monotonically increases with a, for a € (1/2,1). 
This monotonic behavior is an unexpected result that 
strongly contrasts with the uncorrelated case, where the 
gap changes irregularly (discontinuously) with magnetic 
field, as can be inferred from Fig. 1. 

The Hubbard model has already been applied for in- 
vestigations of nanosystems, e^s., molecular wires^ 7 . and 
quantum dot arrays (see Ref. [29| for a detailed discussion 
and estimation of the relevant model parameters). The 



latter case is of particular interest since the size of the el- 
ementary cell can be adjusted in such a way, that the flux 
through the cell can be of the order of the flux quantum 
and, then, the structure of the Hofstadter butterfly may 
be visible. Unfortunately, the present approach is prob- 
ably oversimplified for a quantitative description of elec- 
tronic correlation in the quantum dot arrays, where one 
should account for a large number of states per dot. In 
spite of this, the considered Hubbard Hamiltonian allows 
for at most two electrons per site. Additionally, we have 
restricted our considerations only to the on-site repul- 
sion U, what corresponds to the intradot interaction. In 
a more realistic approach, an extended multiband Hub- 
bard model with intersite interaction should be used, but 
such extensions are presently beyond the reach of the ex- 
act diagonalization techniques. 



IV. SUMMARY 

The structure of the Hofstadter butterfly arises due to 
coexistence of the periodic potential and the perpendicu- 
lar magnetic field in the 2D electron gas. It is well known 
that the presence of Coulomb interaction in low dimen- 
sional systems is responsible for strong modification of 
its electronic properties. Motivated by these facts, we 
have investigated how the butterfly structure is affected 
by the correlations. In contrast to the mean-field results 
we have demonstrated that Coulomb correlations are re- 
sponsible for broadening of the quasiparticle levels. As 
a result some of the Hofstadter bands overlap, and the 
fractal butterfly structure smears out. We have shown 
that for sufficiently strong repulsion the Hubbard gap 
monotonically decreases with the external magnetic field 
for a € (0, 1/2). This result strongly contrasts with the 
irregular behavior of the density of states in the uncorre- 
lated case. Unfortunately, beyond the mean-field level, 
results can be obtained only for relatively small systems. 
We expect, however, that similar changes occur also in 
much larger systems. On the other hand, some of the 
results can directly be applied to the investigations of 
nanosystem, e.g. quantum dot arrays. The finite-size 
effects seriously modify field dependence of its transport 
properties. Similarly to the persistent currents in nano- 
and mesoscopic rings, also the transport currents are dif- 
ferent in systems consisting of odd and even number of 
sites. 
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